Phonon relaxation of subgap levels in superconducting quantum point contacts 
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Superconducting quantum point contacts are known to possess two subgap states per each 
propagating mode. In this note we compute the low-temperature relaxation rate of the upper 
subgap state into the lower one with the emission of an acoustic phonon. If the reflection in the 
contact is small, the relaxation time may become much longer than the characteristic lifetime of a 
bulk quasiparticle. 

In the present paper we address the question of phonon relaxation of subgap levels in superconducting quantum 
point contacts (SQPC) @ §. We use this term for a class of junctions between two superconductors (of the BCS 
type), were only a small number of modes propagates and the scattering center of the contact is much shorter than the 
superconducting coherence length. Under the latter condition, the internal structure of the scatterer is not important, 
but may be described by a scattering matrix for normal electrons Thus, in the class of SQPC we include both 

SNS junctions with a thin normal layer and with quantized propagating modes (^] and mechanically-controllable 
breakj unctions 0. 

Theoretically, SQPCs in different setups are predicted to exhibit various quantum phenomena. The most advanced 
predictions involve quantum-mechanical evolution of a superconducting island connected to external leads by SQPCs 
[§,[)). Other works discuss interference between time evolutions of localized states and/or excitations of localized 
states by an external electromagnetic field |ic|-|l^]. For experimental observation of these predictions, the decay time 
of excited localized states must be sufficiently long. This decay time directly enters the expressions for the thermal 
fluctuations of the Josephson current in SQPC p3|-|T6|. Potential use of SQPC in quantum computing devices also 
crucially depends on their ability to preserve quantum coherence [ p"7| . 

At low temperatures (T <C A, where A is the superconducting gap), the dominant inelastic processes are the 
transitions between localized subgap levels, without participation of continuum spectrum. There are two major 
channels of inelastic relaxation. One is the emission or absorption of a phonon, the other is due to the electromagnetic 
coupling to the external environment. It will depend on the particular experimental setup, which of the two relaxation 
channels dominates. In the present paper we consider only the phonon relaxation in its simplest form, when the 
superconducting phases on the contact terminal are assumed to be rigidly fixed. 

In several works, the time of decay into phonons was estimated as the characteristic time To of a bulk quasiparticle 
p0| |l2[ . The latter lifetime is of order tq ~ Q 2 D / A 3 — the same as the relaxation time of a normal electron with 
energy A above the Fermi level pl| (we set K = 1 throughout the paper). We want to point out that this assumption 
is not justified for decay of localized states. In refs. 0,^9) the relaxation of subgap states in SQPC is associated 
with bulk subgap states appearing due to electron-phonon interaction. This contribution is exponentially small for 
temperatures much lower that the superconducting gap. 

We compute the direct matrix element for decay of localized states into phonons. This direct decay leads to a non- 
vanishing relaxation rate even at zero temperature. Further, we simplify our discussion by setting the temperature 
much lower than the energy of the subgap level. Then the only allowed process is the transition from the upper to 
the lower level with the emission of an acoustic phonon. The extension of our result to include thermal phonons is 
obvious (see eq. (|l6|) below). 

The characteristic energy scale of the subgap levels is A, thus the wavelength of the phonon is of order s/A ~ 
^s/vf *C £, where s and vf are sound and Fermi velocities respectively, £ is the superconducting coherence length. 
The rate and the angular distribution of the emitted phonon may depend on the particular geometry of the contact 
(or, more precisely, on the geometry of the wavefunction of the subgap states). In this paper we discuss the simplest 
setup of a narrow one-dimensional contact. By this we mean that the whole subgap state is localized in a narrow 
strip of width much smaller than the phonon wavelength. Although very idealistic, this assumption is consistent with 
the model of adiabatic constriction H and gives an upper bound for the actual decay rate. 

Each propagating mode in a quasi-one-dimensional contact may be described by the Hamiltonian 



H= dx 
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where and \& are electron operators (L and R subscripts denote left- and right-movers, /3 =|, | is the spin index), 
A (a;) is the superconducting gap with the following x-dependence: 



[It will be convenient for us to use the units with Fermi velocity equal to one throughout the paper] . The scattering 
term H scat t expresses elastic scattering at x — and may be described by a scattering matrix [^j. Diagonalizing the 
Hamiltonian (|l|) gives the subgap state operators 7I and 7] raising energy by 



E(a) = ±A^l-tsm z - (3) 

each (t is the normal transparency of the contact) P,pO|, The two levels below continuum are the ground state |0) 
and the first excited state |1) = t|tJ |0) - The decay of the state |1) to the state |0) with the emission of a phonon 
depends on the density matrix element 

(0|n(x)|l) - (01^(^)^(^)11). (4) 

Since 7! and 7J are linear in electron operators, the matrix element (Q) may be computed by commuting density 
operator with them: 

(0|»(a?)|l) = {[n(a;),7?] = i\b\ K e~ 2 ^sign(x) 7 (5) 

where b is the backscattering amplitude (|6| = VI — t), and k = V~tA \ sin(a/2)| is the inverse length of the subgap 
state. The matrix element is purely imaginary if the relative phases of |0) and |1) are chosen according to § 

(0|^|0) = (1|A| 1) = 0! ( |A|i )isre al. (6) 

This fact is of no importance for the present calculation, but will be used elsewhere in the discussion of the phonon 
emission in the presence of dynamics in a. 

The electron-phonon interaction is described by the deformation potential: 

He- P h = 9 ! d 3 rv(r)*l(r)*p(r), (7) 



tp(r) 



■i{kr-LO k t)\ 



(8) 



where bk are phonon operators normalized by [b^ , b k 



2 = ^C 

1 2e F : 



(in the units with Fermi velocity equal to one), £ is the coupling constant of order one. 
The transition rate is then given by 

1 2 E f cPk 

t- 1 = 2nJ2 \(0,k\H e _ ph \l) 5(<j k - 2E) = ir 3 (— / tt-^K - 2E) (0|n fc |l 



(9) 



(10) 



Here E is the energy of the subgap states given by (||) [so that the energy of the emitted phonon is 2E], Ef is the 
Fermi energy, uik is the phonon dispersion relation, and rik is the three-dimensional density operator. Assume the 
linear isotropic phonon spectrum: 



Wfe = s\k\ 



(11) 
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and the narrow contact limit, where the matrix element of depends only on the component of k parallel to the 
constriction and is given by the Fourier transform of the one-dimensional matrix element (^J). Finally, using k 3> K, 
we arrive to the answer for r _1 : 



T- 1 = 7T 2 C 



E 2 



(ce F ) 



dx 



\n(x)\l) 
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E 2 k 
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Returning to the physical units, we find up to a constant factor of order one 



1 = Vt(i-t) 



(1 



(13) 



where 0£> is the Debye temperature. 

If compared to the characteristic bulk quasiparticle inverse lifetime r^ 1 - A 3 /Q 2 D fL§], the result @ is smaller 
by a factor depending on the backscattering probability 1 — t [|lj. In the case of weak backscattering, this factor 
may contribute up to orders of magnitude to the decay time. This effect is easy to understand: in ideally conducting 
contact (t = 1) the two Andreev states carry opposite momenta equal to the Fermi momentum, and the matrix 
element contains only a rapidly oscillating with momentum 2k F part p3] | : 



\n(x)\l) 



-2k\x\ ^2ikpx 



This oscillating part of the matrix element gives the lower bound for the relaxation rate (p~3h as t — > 1: 



E 3 t 



s 3 e% 
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For realistic values of A, @d, and s F , this relaxation time is unphysically large. The actual relaxation time at t = 1 
will be bounded by other factors such as finite thickness of the interface and non-one-dimensionality of the contact. 
These effects go beyond the simple model of the present paper |£2| . 

Another important feature of the relaxation rate ( |l3| ) is that it vanishes at a — > 0, since in this limit the subgap 
states become delocalized. 

Our assumption of one-dimensionality of the contact also results in overestimating the relaxation rate. If the "tails" 
of the localized states are smeared in the terminals, they give weaker contribution to phonon emission. Unfortunately, 
this effect is highly geometry-dependent, and should be considered separately in each experimental realization. 

As a direct application of the above result, the decay rate entering the fluctuations of the Josephson current in a 
single SQPC OM is given by 



7 = T- 1 (l + 2n B (2E(a)) 



_1 coth 



E{a) 
T ' 



where riB(2E(a)) is the Bose occupation number for the phonons involved in the transition 
temperature rate (|l^). The low- frequency current noise is |I^ , |l(| 



(16) 



L5I , r 1 is the zero- 



s' (a;) = So 



20 



LO 2 + "f 2 
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where Sq is the integral low-frequency noise. 

To summarize, we calculated the rate of direct relaxation of subgap states in SQPC into acoustic phonons at low 
temperature in the simplest one-dimensional geometry. The relaxation rate does not vanish in the T — > limit, but 
it is strongly suppressed in the case of a nearly ballistic contact. 
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# 95-0302, Swiss National Science Foundation collaboration grant # 7SUP J048531, DGA grant # 94-1189, and the 
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The actual lifetime of a bulk quasiparticle diverges at the bottom of the quasiparticle band. The "characteristic" lifetime 
ro corresponds to energies of order A above the bottom of the band. There is no a priori reason for the relaxation rate 
t -1 of the localized states to be of order 

Our results (|l3|), ( |l5| ) disagree with ref. jl5| where the relaxation rate is computed in the same approximation. The 
relaxation rate of ref. jlB) is not suppressed in the ballistic limit t — » 1. 
We neglected the possibility of umklapp scattering. 



4 



